Feigenbaum graphs at the onset of chaos 
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We analyze the properties of networks obtained from the trajectories of unimodal maps at the 
transition to chaos via the horizontal visibility (HV) algorithm. We find that the network degrees 
fluctuate at all scales with amplitude that increases as the size of the network grows, and can be 
described by a spectrum of graph-theoretical generalized Lyapunov exponents. We further define an 
entropy growth rate that describes the amount of information created along paths in network space, 
and find that such entropy growth rate coincides with the spectrum of generalized graph-theoretical 
exponents, constituting a set of Pesin-like identities for the network. 

PACS numbers: 05.45.Ac, 05.45.Tp, 89.75.Hc 



I. INTRODUCTION 



Pesin's theorem [l| prescribes the equality of the 
Kolmogorov- Sinai (KS) entropy has with the sum of 
the positive Lyapunov exponents Ai > of a dynam- 
ical system, i.e. h^s — Y)j Aj. The former is a rate 
of entropy growth that is a metric invariant of the dy- 
namical system, while the latter is the total asymptotic 
expansion rate present in the chaotic dynamics. This re- 
lation provides a deep connection between equilibrium 
statistical mechanics and chaos. In the limiting case of 
vanishing Lyapunov exponent, as in the onset of chaos 
in one-dimensional nonlinear maps, Pesin's theorem is 
still valid but there are important underlying circum- 
stances that are not expressed by the trivial identity 
fiKS = A = 0. At variance with the chaotic region, at 
the transition to chaos phase space is no longer visited in 
an ergodic way and trajectories within the attractor show 
self-similar temporal structures, they preserve memory of 
their previous locations and do not have the mixing prop- 
erty of chaotic trajectories For chaotic and periodic 
attractors the sensitivity to initial conditions converges 
to an exponential function for large iteration times giv- 
ing rise to the familiar positive or negative Lyapunov 
exponents associated with them. At the transition to 
chaos exponential separation or merging of trajectories 
no longer occurs, but notably the sensitivity to initial 
conditions does not converge to any single-valued func- 
tion and, on the contrary, displays fluctuations that grow 
indefinitely with time. For initial positions on the attrac- 
tor the sensitivity develops a universal self-similar tem- 
poral structure and its envelope grows with iteration time 
t as a power law @, Q . For unimodal maps it has been 
shown [2|, Q that this rich borderline condition accepts a 
description via a spectrum of generalized Lyapunov ex- 
ponents that matches a spectrum of generalized entropy 
growth rates obtained from a scalar deformation of the 
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ordinary entrow functional, the so-called Tsallis entropy 
expression 0, |f|. The entropy rate also differs from the 
KS entropy in that it is local in iteration time. That is, 
the entropy rate at time t is determined only by the po- 
sitions of trajectories at this time while the KS entropy 
considers all previous sets of positions. Here we present 
evidence that this general scenario of unimodal maps at 
the onset of chaos is captured by a special type of com- 
plex network. 

Very recently [f| 0] , the horizontal visibility (HV) algo- 
rithm that transforms time series into networks has 
offered a view of chaos and its genesis in low-dimensional 
maps from an unusual perspective favorable for the dis- 
covery of novel features and new understanding. We fo- 
cus here on networks generated by unimodal maps at 
their period-doubling accumulation points and charac- 
terize the fluctuations in connectivity as the network size 
grows. We show that the expansion of connectivity fluc- 
tuations admits the definition of a graph-theoretical Lya- 
punov exponent. Furthermore, the entropic functional 
that quantifies the amount of information generated by 
the expansion rate of trajectories in the original map 
appears to translate, in the light of the scaling proper- 
ties of the resulting network, into a generalized entropy 
that surprisingly coincides with the spectrum of gener- 
alized graph-theoretical Lyapunov exponents. This sug- 
gests that Pesin-like identities valid at the onset of chaos 
could be found in complex networks that possess certain 
scaling properties. 

The rest of the Letter is as follows: We first recall the HV 
algorithm that converts a time series into a network and 
focus on the so-called Feigenbaum graphs [1, 0] as the 
subfamily of HV graphs generated by iterated nonlin- 
ear one-dimensional maps. We then expose the univer- 
sal scale-invariant structure of the Feigenbaum graphs 
that arise at the period-doubling accumulation points. 
We subsequently show how to define a graph-theoretical 
Lyapunov exponent in this context. In agreement with 
the known dynamics of unimodal maps at the transition 
to chaos the graph-theoretical Lyapunov exponent van- 
ishes, and proceed to define generalized exponents that 
take into account the subexponential expansion of con- 
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FIG. 1: Feigenbaum graphs associated with periodic series of 
increasing period 2" undergoing a period-doubling cascade. 
The resulting patterns follow from the universal order with 
which an orbit visits the positions of the attractor. The 
Feigenbaum graph associated with the time series generated 
at the onset of chaos (n — > oo) is the result of an infinite 
application of the inflationary process by which a graph at 
period 2 n+1 is generated out of a graph at period 2" 



nectivity fluctuations. We finally show that the Feigen- 
baum graph that represents the onset of chaos admits a 
spectrum of generalized graph-theoretical exponents that 
coincide with a spectrum of deformed entropies. 



II. THE FEIGENBAUM GRAPH AT THE 
ONSET OF CHAOS 

The horizontal visibility (HV) algorithm is a general 
method to convert time series data into a graph [1, HI 
and is concisely stated as follows: assign a node i to 

each datum Xi of the time series {xt\i = \ n of N real 

data, and then connect any pair of nodes i, j if their 
associated data fulfill the criterion Xi, Xj > x n for all n 
such that i < n < j. The capability of the method as 
well as other variants to transfer properties of different 
types of time series into their resultant graphs has been 
demonstrated in recent works [lfj, . When the series 
under study are the trajectories within the attractors 
generated by unimodal or circle maps the application 
of the HV algorithm yield subfamilies of visibility 
graphs, named Feigenbaum and Quasiperiodic graphs 
respectively, that render the known low-dimensional 
routes to chaos in a new setting @, 0, [H[ • For illustrative 
purposes, in Fig. [TJ we show a hierarchy of Feigenbaum 
graphs obtained along the period-doubling bifurcation 
cascade of unimodal maps. At the accumulation point 
of the cascade the trajectories within its non-chaotic 
multifractal attractor become aperiodic, and in analogy 
with the original Feigenbaum treatment [13], the associ- 
ated Feigenbaum graphs evidence scaling properties that 
can be exploited by an appropriate graph-theoretical 
Renormalization Group transformation [6j, [7]. Because 
the order of visits of positions of periodic attractors or 



of bands of chaotic attractors in unimodal maps are 
universal, a relevant consequence of the HV criterion 
is that the resulting Feigenbaum graph at the onset of 
chaos is the same for every unimodal map. That is, it is 
independent of the shape and nonlinearity of the map 
0, Q. This permits us to concentrate our study on a 
specific unimodal map, that for simplicity will be the 
logistic map, and claim generality of the results. 

Consider the logistic map xt+i = f(%t) = ~~ x t), 

0<a:<l,0</i<4, at the period-doubling ac- 
cumulation point //oo = 3.5699456... and generate a 
trajectory {xt}t=o.i,2.... within the attractor. Since 
the position of the maximum of the map x ma x = 1/2 
belongs to the attractor we can choose xq = /(1/2) as 
the initial condition. (The position of the maximum 
of the unimodal map belongs to all the superstablc 
periodic attractors or supercycles and the accumulation 
point of these attractors, the attractor at the transition 
to chaos, contains this position [HI])- A sample of the 
corresponding time series is shown in the top panel of 
figure [2j whereas in the left panel of the same figure 
we represent in logarithmic scales the positions of the 
rescaled variable z t — /(1/2) — Xt- A similar rescaling 
was proposed in previous works to put in evidence 
the striped intertwined self-similar structure of the 
trajectory positions as they are visited sequentially, and 
reflects the multifractal structure of the attractor 0, H| . 
In the right panel of the same figure we plot, in log-log 
scales, a representation of the relevant variable exp k(N), 
where k(N) is the degree of node N in the Feigenbaum 
graph associated to the original time series x(t) (that is, 
N = t). Notice that the distinctive band pattern of the 
attractor is recovered, although in a simplified manner 
where the fine structure is replaced by single lines 
of constant degree. The HV algorithm transforms the 
multifractal attractor into a discrete set of connectivities, 
whose evolution mimics the intertwined fluctuations 
of the map. It is convenient to write the node index 
N in the form N = m2 J , where j = 0,1,2,... and 
m = 1, 3, 5, so that running over the index j with m 
fixed selects subsequences of data or nodes placed along 
lines with fixed slope. For illustrative purposes, the 
order of visits to either the associated data or nodes in 
some of these subsequences are explicitly shown in both 
panels of figure O as well as in the upper panel. Notice 
that along these subsequences, expfc(TV) increase as a 
power law (i.e. subexponentially) , and are associated to 
data with increasing larger values (upper panel). This is 
reminiscent of the data subsequences highlighted in the 
left panel, whose values z(t) = /(1/2) — x(t) also show a 
monotonic power law decay. 

From the direct observation of figure [H it is clear that 
all trajectory positions within each band in the left 
panel (attractor) become nodes with the same degree 
in the right panel (network). It is also worth noting at 
this point that there exists a degeneracy in the positions 
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FIG. 2: Up: Series x(t) as a function of time t for the first 10 6 data generated from a logistic map at the period-doubling 
accumulation point (only the first 33 data are shown). The data highlighted in red are associated with specific subsequences 
of nodes (see the text). Left: Log-log plot of the rescaled variable z(t) = /(1/2) — x(t) as a function of t, for the same series as 
the upper panel. This rescaling is performed to reflect the multifractal structure of the attractor 0, 0] • The order of visits to 
some specific data subsequences is highlighted. Right: Log-log plot of exp k(N) as a function of the node N of the Feigenbaum 
graph generated from the same time series as for the upper panel, where N — t. The distinctive band pattern of the attractor 
is recovered, although in a simplified manner where the fine structure is replaced by single lines of constant degree. The order 
of visits to some specific node subsequences is highlighted (see the text). 



of trajectories that are initiated off, but close to, the 
attractor positions. In particular, if an ensemble of 
uniformly-distributed initial positions is placed in a 
small interval around x = 1/2 their trajectories expand 
while remaining uniformly-distributed at later iteration 
times 2, 3] (see discussion below) and the HV algorithm 
assigns to all of them the same Feigenbaum graph. This 
degeneracy accounts for the universal feature that only 
a single Feigenbaum graph represents the transition to 
chaos for all unimodal maps. It also determines the way 
in which fluctuations in trajectory separations translate 
into their respective degree fluctuations: see figure [3j 
where we plot in logarithmic scales the distance of 
initially nearby trajectories for specific values of a time 
subsequence, and compare the expansion of this distance 
with the increment of the degree of the associated 
Feigenbaum graph along the same node subsequence, as 
highlighted in the right panel of figure [2] The properties 



of trajectory separation is thereby inherited by the graph 
and can be quantified through the increasing sequences 
of values for exp k(N). 

Let us focus now in the Feigenbaum graph. The 
relabeling N = mV generates a one-to-one tiling of the 
natural numbers. The aforementioned constant degree 
and fixed-slope lines in the right panel of figure [5] are, re- 
spectively, k(N = m2^) = 2j + 2, m = 1, 3, 5, j fixed, 
and k(N = m2 ] ) = 2j + 2, j = 0, 1,2, m fixed (for 
example, the first path highlighted in red corresponds to 
m = l,j = 0,1,2,...). The scaling property associated 
with the former case is that all horizontal lines with 
constant degree k = 2 j + 2, j fixed, can be overlapped 
into the bottom line k = 2, j = 0, via consecutive 
translations each consisting of a shift of 2 in k and a 
shift of log 2 in log TV, that is, k(N/2) = k(N) - 2, 
or alternatively, expk(N/2) = exp(— 2) • expfe(A^). 
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FIG. 3: Log-log plot of the distance between two nearby 
trajectories d(t) = \x(t) — x'(t)\ close to x = 1/2, where 
d(0) = 10 -6 , measured at a precise time sequence t = 2^,j = 
1,2,3,... (this path is also highlighted in the left panel of 
figure [2j . Along this path trajectories evidence a subexpo- 
nential expansion characterized by a power law separation. 
Note that both trajectories generate the same Feigenbaum 
graph, whose degree increases along this concrete path in the 
same vein as the trajectory separation: the network inherits 
the properties of trajectory expansion in the hierarchy of de- 
grees (see right panel of figure [2} . A similar behavior holds 
for trajectory separation along other paths. 



This scaling property has a parallel for the trajectory 
within the Feigenbaum attractor where the shift of 2 
in k is replaced by a shift of 2 log a, where a is the 
absolute value of the Feigenbaum constant [H, The 
scaling property associated with the latter case is seen 
via the collapse of all fixed-slope lines of data into a 
single sequence of values aligned along the first m = 1 
subsequence fc(2 3 ) = 2j + 2, j = 0,1,2,..., when the 
numbers of nodes for each of the other subsequences are 
rescaled consecutively by N/m, m — 3, 5, 7, .... 



III. FLUCTUATING DYNAMICS AND 
GRAPH-THEORETICAL LYAPUNOV 
EXPONENTS 

We recall that the standard Lyapunov exponent A ac- 
counts for the degree of exponential separation of nearby 
trajectories, such that two trajectories whose initial sepa- 
ration is d(0) evolve with time distancing from each other 
as d(t) ~ d(0) exp(At), t » 1, where d(t) = \x t -x' t \. The 
sensitivity to initial conditions is = d(t)/d(0) and 
the general definition of the Lyapunov exponent is 



Deterministic fluctuations about A > are permanently 
put out when t — > oo, but, as mentioned, at the period- 
doubling accumulation points of unimodal maps A = 
and trajectories display successive subexponential sepa- 
ration and convergence for all t. Two possible extensions 
of the standard Lyapunov exponent have been used when 
A = while the amplitude of the fluctuations grows as 
a power law, or, equivalently, expansion and contraction 
rates of trajectories are logarithmic in time: Mori et al 
13] suggested the expression 
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log* 



log£(f). 



(2) 



Although this expression captures subexponential fluctu- 
ation rates, this quantity is not well defined if we want to 
measure the magnitude of rates per unit time, in which 
case we need a quantity that grows linearly in time. This 
is important if we want to relate the growth of fluctua- 
tions to some kind of entropic rates that derive from an 
extensive entropy expression. In this respect a better- 
suited alternative for capturing subexponential fluctu- 
ations is to deform the logarithm in equation Q] by an 
amount that recovers the linear growth present for A > 0. 
This is 



(3) 



with \og q x = [x q —l)/(l — q) (logx is restored in the 
limit q — > 1) where the extent of deformation q is such 
that while the expansion (ft) is subexponential, log_£(t) 
grows linearly with t 0, Q ■ Notice that we have elimi- 
nated the t — > oo limit in the definition of the generalized 
Lyapunov exponents since fluctuations are present for all 
values of t and the objective is to characterize them via 
a spectrum of generalized exponents @, d, [l]| . 

We define now a connectivity expansion rate for the 
Feigenbaum graph under study. Since the graph is con- 
nected by construction (all nodes have degree k > 2), 
the uncertainty is only associated with a rescaled de- 
gree k+ = k — 2. To keep notation simple we make 
use of this variable and drop the subindex + from now 
on. The formal network analog of the sensitivity to ini- 
tial conditions has as a natural definition £(N) = 
exp k(N) j exp fc(0) = expk(N), as suggested from fig- 
ures [21 and where we are implicitly assuming that 
the expansion is always compared with the minimal one 
exp/c(0) = 1 provided by nodes at positions to2° (i.e. 
with k = 0). Accordingly, the standard network Lya- 
punov exponent is defined as 



A= lim -jzlogtfJV). 



(4) 



A= lim -log£(t). 

t->oo t 



(1) 



The value of k(N) oscillates with N (see figure [5]) but its 
bounds grow slower than TV, as logiV, and therefore in 
network context A = 0, in parallel to the ordinary Lya- 
punov exponent at the onset of chaos. The logarithmic 
growth of the bounds of log£(AT) = k(N) is readily seen 
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by writing k(N = m2 J ) = 2j as 

k(N) = -^- log (- 
log 2 \m 



(5) 



A first approach to study the subexponential fluctuations 
is to proceed a la Mori and define the following graph- 
theoretical generalized exponent 



A 



1 



M 



log TV 



log£(JV) - 



k(N) 
logiV' 



(6) 



To visualize the network growth paths A — m2 J , where 
m = 1,3,5,... is fixed to a constant value and j = 
1, 2, 3, the expression above for Am is written as 



Am — 



2j 



2j 



log(m2J ) log m + j log 2 



(7) 



A constant Am = 2/ log 2 for all j is obtained only when 
m = 1, otherwise the same value is reached for all other 
m when j —t oo. In general, a spectrum of Mori- like 
generalized exponents is obtained by considering all 
paths (m,j) such that A — > oo. 

As indicated, the preceding generalization of the Lya- 
punov exponent is not time extensive and it is therefore 
not useful if we ultimately wish to relate the network 
link fluctuations to an entropy growth rate. With this 
purpose in mind we deform the ordinary logarithm in 
log£(A) = k(N) into log ? £(iV) by an amount q > 1 such 
that log^ £,(N) depends linearly in N, and define the asso- 
ciated generalized graph-theoretical Lyapunov exponent 
as 



(8) 



where AA = N — m is the lapse time between an ini- 
tial condition m2° where m is fixed (see figure [5) , and 
position A. One obtains 



~~ — i — o' 
m log 2 



(9) 



with q = l— log 2/2 and j > 0. Eq. (8) can be corrobo- 
rated via use of £(A) = exp(2j) together with the iden- 
tifications of the other quantities in it (note that j ' = 
implies AA = for which £(A) = exp(0) = 1 and X q 
is trivially undefined). According to Eq. (9) a spec- 
trum of exponents is spanned by running over the values 
of m, in parallel with the spectrum of generalized ex- 
ponents previously found at the transition to chaos in 
unimodal maps [3j, where the value of the parameter, 
q = 1 — log 2/2 in the Feigenbaum graph is to be com- 
pared with q = 1 — log 2 / log a (where a is a Feigenbaum 
constant) in the unimodal map for trajectories originat- 
ing at the most compact region of the multifractal at- 
tractor that are seen to expand at prescribed times when 
the least compact region of the multifractal is visited Q . 



IV. ENTROPIC FUNCTIONALS AND 
PESIN-LIKE IDENTITIES 

To complete our arguments we summon up the per- 
sistency property of trajectory distributions of unimodal 
maps at the period-doubling onset of chaos. That is, for a 
small interval of length l\ with M uniformly-distributed 
initial conditions around the extremum of a unimodal 
map {e.g. x = 1/2 for the logistic map we use), all tra- 
jectories behave similarly, remain uniformly-distributed 
at later times, follow the concerted pattern shown in the 
left panel of Fig. 2. (see 0, Q for details) and expand 
subexponentially at prescribed time paths (see figure [3]). 
Studies of entropy growth associated with an initial dis- 
tribution of positions with iteration time t of several 
chaotic maps [lj| have established that a linear growth 
occurs during an intermediate stage in the evolution of 
the entropy, after an initial transient dependent on the 
initial distribution and before an asymptotic approach 
to a constant equilibrium value. At the period-doubling 
transition to chaos it was found that (i) there is 

no initial transient if the initial distribution is uniform 
and defined around a small interval of an attractor po- 
sition, and (ii) the distribution remains uniform for an 
extended period of time due to the subexponential dy- 
namics. We denote this distribution by n(t) = 1/W(t) 
where W(l) = l\/N is the number of cells that tile the 
initial interval l\, and It is the total length of the interval 
that contain the trajectories at time t. Since iteration- 
time dependence transforms into node- value dependence 
we inquire about how 7r scales along the set of nodes 
A = V with k = 2j links, j = 0, 1, 2,. . . , that is, while 
7r is defined in the map its properties are reflected in the 
network. After t = 2 J iterations, the loss of information 
associated with the fact that the same network is gener- 
ated by all of the M initial conditions is represented in the 
network context by the expansion of the initial interval 
li into l t = exp fc(2 J ) = li exp(2j) after t = 2P. Therefore 
we find that at time t = 2 J the network scaling properties 
prescribe a larger number of cells to tile the new interval 
W{2^) = W(l) exp fc(2- J ). After normalization, we find 
that the j dependence of 7r reads 



7r(2 J ') = W; 1 = exp(-2j). 



(10) 



That is, the uniform distributions -k for the consecutive 
node-connectivity pairs (2 J , 2j) and (2 J+1 , 2(j + 1)) scale 
with the same factors noticed above when discussing the 
right panel of figure 2. Namely, when m ~ 1 these pairs 
can be made equal by a shift of log 2 in log A and a 
shift of 2 in k. By extension of the argument the same 
expression holds for all other values of m. Since 



Wj = exp(2j) 



m J 



(11) 



the ordinary entropy associated with n grows logarithmi- 
cally with the number of nodes A, S\ [7r(A)] = log Wj ~ 



6 



log N. However, the g-deformed entropy 



S q [ir(j)}=log q W j = - 



1 



Q 
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(12) 



where the amount of deformation q of the logarithm has 
the same value as before, grows linearly with N, as Wj 
can be rewritten as 



Wj — expJAqAiV], 



(13) 



with q = 1 - log 2/2 and A, = 2/ (m log 2). Therefore, if 
we define the entropy growth rate 



h q [tt(N)} 



1 



N -m 



S q [tt(N)} 



we obtain 



h q [tt(N)} = X q , 



(14) 



(15) 



a Pesin-like identity at the onset of chaos (effectively 
one identity for each subsequence of node numbers given 
each by a value of m — 1, 3, 5, ...). 



V. CONCLUSIONS 

In conclusion, the transcription into a network of a spe- 
cial class of time series, the trajectories associated with 
the attractor at the period-doubling transition to chaos 
of unimodal maps, via the HV algorithm has proved to be 
a valuable enterprise @ as it has led to the uncovering of 
a new property related to the Pesin identity in nonlinear 
dynamics. The HV method leads to a self-similar net- 
work with a structure illustrated by the related networks 
of periods 2", n = 0,...,5, shown in Fig. [T] Under the 
HV algorithm many nearby trajectory positions lead to 
the same network node and degree, all positions within 
one band in the left panel of Fig. [2] lead to the same line 



of constant degree in the right panel of the figure. Only 
when trajectory positions cross a gap between bands in 
the left panel the corresponding node increases its de- 
gree by two new links. Also trajectories off the attractor 
but close to it transform into the same network struc- 
ture. The degrees of the nodes span all even numbers 
k = 2j, j = 0,1,2,..., and we have studied how these 
fluctuate as the number of nodes increases. The fluc- 
tuations of the degree capture the core behavior of the 
fluctuations of the sensitivity to initial conditions at the 
transition to chaos and they are universal for all unimodal 
maps. The graph-theoretical analogue of the sensitivity 
was identified as exp(fc) while the amplitude of the varia- 
tions of k grows logarithmically with the number of nodes 
N. These deterministic fluctuations are described by a 
discrete spectrum of generalized graph-theoretical Lya- 
punov exponents that appear to relate to an equivalent 
spectrum of generalized entropy growths, yielding a set 
of Pesin-like identities. The definitions of these quan- 
tities involve a deformation of the ordinary logarithmic 
function that ensures their linear growth with the num- 
ber of nodes. Therefore the entropy expression involved 
is extensive and of the Tsallis type with a precisely de- 
fined index q. A salient feature of the application of the 
HV algorithm to dynamical systems time series is direct 
access to the degree distribution and therefore to the en- 
tropy associated with it. Therefore this seems to be a 
good tool to discover, or, inversely, to construct, net- 
work entropy properties as studied here. The capability 
of the methodology employed in this study to reveal the 
occurrence of structural elements of this nature in real 
networked systems that grow in time [l6l . [l7j appears to 
be viable. 
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